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We study transverse and longitudinal electron heating effects on the target stability and the ion 
spectra in the radiation pressure dominated regime of ion acceleration by means of multi dimensional 
particle-in-cell (PIC) simulations. Efficient ion acceleration occurs when the longitudinal electron 
temperature is kept as low as possible. However, tailoring of the transverse electron temperature is 
required in view of suppressing the transverse instability, which can keep the target structure intact 
for longer duration during the acceleration stage. We suggest using the surface erosion of the target 
to increase the transverse temperature, which improves both the final peak energy and the spectral 
quality of the ions in comparison with a normal flat target. 
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In recent years, ion acceleration from thin foil targets 
has emerged as one of keyareas of research in the field of 
laser plasma interaction [l| . The laser foil target interac- 
tion can produce energetic ions with energies as high as 
56MeV per nucleon. Target engineering is crucial to get 
the mono-energetic ion beams [2| . These ion beams have 
extremely short duration (~fs), highly collimated and are 
relatively easy to produce, which makes them suitable for 
many applications, such as proton imaging (sj, ion ther- 
apy Q , ion beam ignition of laser fusion targets [H and 
so on. It was recently suggested to utilize the circularly 
polarized (CP) laser pulses for very high energy ion accel- 
eration [6]. In the case of CP pulse, the electron heating 
is dramatically reduced, thus ion acceleration doesn't oc- 
cur due to the target normal sheath acceleration (TNSA), 
instead radiation pressure acceleration (RPA) dominates. 
Though energy scalings predicted by ID theory of RPA 
regime of acceleration have been reproduced very well in 
numerous one dimensional particle-in-cell (PIC) simula- 
tions, the accelerated target is far less stable in a real 
3D geometry compare to ID geometry. For instance, in 
the case of a Gaussian pulse interacting with a flat target, 
target deformation occurs, which leads to the broadening 
of the final energy spectrum, and reduces the energy con- 
version efficiency; though it can be overcome by the use 
of shaped targets [7]. Another major bottleneck in the 
ion acceleration is the excitation of the Ray leigh- Taylor 
like instabilities, which leads to the breaking of the tar- 
get [8j,5]. Thus, efficacy of the RPA scheme is limited by 
the onset of the Ray leigh- Taylor like instability; which 
also appears to be a cause of great concern in ion beam 
driven fast ignition. Prevention of this instability is an 
important problem which needs to be addressed quickly. 

As described above, the RPA mechanism dominates 
when the electron heating is suppressed; this condition 
can be fairly met if one resorts to the use of the CP laser 
pulses. However, the onset of the Ray leigh- Taylor like 
instability, also known as transverse instability, remain 
a concern for the CP laser pulses. With regards to elec- 
tron temperature, there is tradeoff as large perpendicular 
temperature might suppress the Ray leigh- Taylor like in- 



stability while longitudinal electron temperature can lead 
to the broadening of the energy spectrum. In this paper, 
we explore such possibilities and study the heating ef- 
fects by differentiating the transverse and longitudinal 
heating of the electrons. We find that for efficient ion ac- 
celeration, longitudinal heating of the electrons should be 
kept as low as possible, however, proper transverse heat- 
ing can suppress the transverse instability and keep the 
target structure uniform for longer durations. Since the 
CP laser pulse can reduce both transverse electron tem- 
perature Te±_ and longitudinal electron temperature Tgy, 
we suggest using target surface erosion to increase Tgy. 
We demonstrate, through multi-dimensional PIC simula- 
tions, that this active controllable method can keep the 
target structure uniform for longer time compared with a 
normal flat target. We believe that our differential treat- 
ment of the transverse and longitudinal electron heating 
not only benefits the ion acceleration but could also be 
useful in the fast ignition scheme. 

We first show the results of 2D-PIC simulations per- 
formed by using the VLPL-code [10]. The size of the 
simulation box is 15Ao(x) x2Ao(^) with Aq representing 
the laser wavelength, which corresponds to a grid of 
3000 (x)x 200 (?/). The time step of the simulation is 
0.003To, here Tq = 3.33 fs is the laser period. The foil tar- 
get consists of two species, electrons and protons. They 
are initially located in the region 2Ao < x < 2.2Ao with 
the density of n = 160nc, where Uc = cj^me/47re^ is the 
critical density for the laser pulse with the frequency u; 
which is ric = 1.1 X 10^^ /cm^ for l/j,m wavelength laser 
pulse. We use 150 particles per cell to run the simula- 
tions. To exclude other effects arising due to the shape 
of the laser pulse, we here only take the plane wave laser 
pulse. The normalized amplitude of the laser electric field 
is ao = eEo/rrieUJC = 100. The laser pulse has a trape- 
zoidal temporal intensity profile (linear growth - plateau - 
linear decrease) , with IAq/c — 8Ao/c— IAq/c. Att = the 
laser pulse enters the simulation box from the left bound- 
ary. Two kinds of targets are used: a normal flat target 
with initially zero temperature, and a target with sur- 
face erosion and initially zero temperature. The second 
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Figure 1: (color online) Spatial density distribution of elec- 
trons (left column) and protons (right column) at t = 20To. 
(a,b) a flat target with initially zero temperature; (c,d) a mod- 
ulation target (Idepth — 0.05Ao,Ay = 0.4Ao) with initial zero 
temperature; (e,f) a flat target with initial transverse temper- 
ature of Te = 6.3Mey 
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Figure 2: (color online) (a) Proton spectrum in the cases 
of a flat target with thickness of 0.2 Ao (black solid line), 
0.15 Ao (green solid line) and a target with modulation (red 
solid line) at t = 20To. (b) Longitudinal (Tex , black line) 
and transverse (Teyz, red/grey line) electron temperatures at 
t = 5 To (solid line) and t = 20 To (dashed line) for a flat tar- 
get and initial zero temperature, (c) Longitudinal (Tex , black 
line) and transverse (Teyz, red/grey line) electron tempera- 
tures at t = 5To (solid line) and t = 20To (dashed line) for 
a target with modulation of Xy = 0.4Ao (black solid line) and 
initial zero temperature, (d) Proton spectrum, at t = 20 To, in 
the cases of initially longitudinal (black solid line) and trans- 
verse temperature (red solid line). 



target has a ripple on the left surface, whose boundary 
is defined as x{y) = xq — ldepthx[s'^^{'^'7ry / Xy) + l]/2 and 
Idepth = 0.05Ao,Ay = 0.4Ao. 

Fig. [T] shows the spatial target density distribution at 
t = 2OT0. As we can see, the normal flat target [sub- 
plots (a) and (b)] is almost completely dissociated after 
20 laser periods. The target is longitudinally expanded 
and transversely agglomerated. The flat target with an 
initial surface ripple [subplots (c) and (d)] has a better 
profile. Both the electrons and protons have sharp longi- 
tudinal boundaries and uniform transverse distribution. 
The main target has been longitudinally compressed to 
a layer of O.OTAo thickness. Its central position is 0.3Ao 
forward than the normal flat target. The proton spec- 
trums are shown in Fig. [2j The spectrum of the normal 
flat target is shown in Fig. [2fa) with the black solid line, 
while red line depicts the spectrum of the eroded target. 
It's clear that both the peak energy and width of the 
spectrum have been improved by using the surface ero- 
sion of the target. To exclude the target thickness effect, 
we also show the spectrum of a cold flat target with a 
thickness of 0.15 A with the green solid line. It also shows 
a broaden spectrum whose peak energy is smaller than 
the eroded target, which has a larger mass. Thus, it is 
clear that the eroded target can keep the acceleration 
structure intact for longer duration, thereby increasing 
the peak energy and reducing the width of the spectrum 



of the accelerated protons. 

In order to ascertain the effect of electron temperature 
on the proton spectrum, we plot the target electron mo- 
mentum distribution at different times for both targets. 
The electron temperature of the flat target and the sur- 
face eroded target at t = 5To (solid line) and at t = 20To 
(dashed line) are shown in Fig.[2](b) and (c), respectively. 
During the early stage of acceleration, the electrons in the 
flat target acquire lower longitudinal temperature {Tex) 
as well as transverse temperature {Teyz) compared with 
the surface eroded target. This implies that the target 
has been well compressed and the onset of instability has 
not occurred, resulting in target being pushed uniformly 
in the forward direction. The electrons in the surface 
eroded target attain higher temperatures perhaps due to 
higher absorption of the laser energy-caused by the sur- 
face ripples- into the target. However, the longitudinal 
temperature is not high enough to break the target as 
the electrons' energies are far less than the ponderomo- 
tive potential level ( a^^GeV), and most of the electrons 
are within the target area. However the scenario changes 
and at later times 20To, we see both of the transverse 
and longitudinal electron temperatures are higher in the 
flat target than those of the surface eroded target. The 
higher electron temperatures in the flat target case are 
due to the breaking and subsequent heating of the target 
by the laser pulse. However, the transverse electron tem- 
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perature of the surface eroded target hasn't witnessed 
higher growth. Most of the electrons in the target are 
around the proton layer to form the uniform high den- 
sity plasma layer to suppress the laser penetration as also 
depicted in Fig. [T] [subplot (c)]. One may infer during 
the early stage of acceleration (T = 5To), the higher elec- 
tron temperature in the surface eroded target appears to 
have provided a stabilizing influence on the transverse in- 
stabilities of the target, thus preventing the target from 
breaking. 

In order to further isolate the temperature effects on 
the target structure during the acceleration stage, we 
perform another simulation by taking a fiat target with 
an initial transverse temperature of 6.3MeV. The tar- 
get electron and proton density distributions are shown 
in Fig. [TJe) and (f). Again, we see a well accelerated 
target foil. The center of the target is little backward 
compared with the erosion target, which is due to the 
heavier mass of the target in this case. The proton spec- 
trum at t = 20To is shown in Fig. [2fd) with the red 
sohd line. It shows a narrow spectrum compared with 
the normal cold flat target case. We also check the ef- 
fect of the longitudinal temperature by using a fiat tar- 
get with an initial longitudinal temperature of 6.3MeV. 
This temperature is not high enough to affect the RPA 
acceleration. However the final spectrum [see the black 
solid line in Fig. Efd)] is even worse then the one in the 
normal fiat target case. This means the longitudinal tem- 
perature cannot suppress the transverse instability and 
on the contrary it increases it. 

From analytical view, to give a quantitative expla- 
nation one should use multi dimensional relativistic 
warm plasma theory with anisotropic relativistic tem- 
perature. For the moment, such kind of theory is 
still not available. Here we give a qualitative descrip- 
tion by reducing the problem to two dimensional ge- 
ometry and seeing the balance between the transverse 
temperature pressure and the perturbed light pressure. 
In the following we use normalized variables, such as: 
x^kpX^t^ujpt^n^n/nc^p^p/meC. kp^uOp^Uc are the 
normal plasma wave number, frequency and critical den- 
sity, respectively. We assume a circularly polarized laser 
pulse interacting with an initial fiat target whose trans- 
verse relativistic temperature is T^y. The density per- 
turbation is in the transverse direction (in y direction) 
n = no + (5n, and 5n = ni sin(27r7//A^)/(x) is the sec- 
ond density perturbation term. The initial uniform elec- 
tron density is n^. For a short time {t « 1/ujpi), ions 
are assumed to be immobile. We will see the transverse 
temperature pressure can weaken the density perturba- 
tion due to light pressure. Because of two dimensional 
condition (x-y plane), we assume all the physical vari- 
ables satisfy dz = 0. For a relativistic warm collision 
plasma, one can use Vlasov equation to derive the fluid 
equations. Here we use the fundamental results given 
by Schroeder et a/. [ill] and use further assumption to 



simplify our model. Although our temperature is rel- 
ativistic ( 6MeV) we still adopt Schroeder's model to 
give a qualitative description. In further, in our case 
jth << \/l + ^0 (Here ^yth corresponds to the temper- 
ature relativistic factor in Schroeder's paper and ao is 
the laser intensity we used here.), we neglect the thermal 
correction to the relativistic Lorentz factor as shown in 
formula (30) in Ref.fl]]]. The second term of relativis- 
tic temperature pressure in formula (30) has also been 
neglected due to time averaging over the laser cycle. 

Then from the anisotropic warm plasma fluid equation, 
one can get the first and second order of the laser cycle 
averaged electron fluid equations as following. 

The continuous equation gives: 
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We have already assumed vqx = and dvoy/dy = 0, 
dno/dy = here, then we can get dno/dt = 0, which 
shows the first order of balance between the light pres- 
sure and the static electric field resulting from the charge 
separation. However, the second order of density Sn is 
the density perturbation due to the perturbed light pres- 
sure. The Poison equation gives: 



no - ni 
Sn 



For the motion of electron fluid, we have: 

dp p da p ^ V ■ {j'^n < UtiUtj >) 

— + --Vp = — + V^--x VxaH ^ — - 

at J ot J n 

(1) 

Here uu is the normalized electrons' thermal momentum 
in the direction, the bracket means taking the average 
over its distribution, n is the proper density of the fluid el- 
ement. Considering the first order of motion (which gen- 
erates the main electron density perturbation: no — n^), 
since we have assumed dy = 0^ for the longitudinal di- 
rection we have: [V • (7^n < uuUtj >)]\x = ^^(To^o < 
utxUtx >) =< > dxijo^o) ^lo^odx < ulx > /2, 
where the underline means the derivation only works 
on that variable. For the transverse direction we have 
[V • (7^n < utiUtj >)]\y = dxijo^o < utxUty >) = 0, 
since < UtxUty >=< Utx >< Uty >= when the ther- 
mal momenta in different directions are not correlated. 
For our interests, < ul^ > <C < uly > (equally to say 
Tex<^Tey)^ thus we assume the longitudinal temperature 
effect could also be neglected, so in the first order of 
electron motion there is no temperature effect. The first 
order of electron motion equation is as following: 
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Figure 3: (color online) Schematic profile of laser foil interac- 
tion. The red arrow shows the laser pressure and perturbed 
light pressure, the transparent blue arrow shows the thermal 
pressure. 



In further we have: 

^(Po - ^o) 

dt 



V(^o - 7o) 



And po± ao±, pox ^ox = 0, 70 = \/l + a§. 
The second order of electron motion equation is: 



d{Sp — 5a) 
di 



- J7) + — X V X {5p- 5a) 
7o 
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no 

For the last term, again we have: {V • [(7o^n + 
27ono^7) < UtiUtj >]}\x = dx[{ll5n + 27ono^7) < 
UtxUtx >]^dy[{^l5n^ 2^071^5^) < UtyUtx >] ^ since < 
utx >— and < u^^ >c^ 0, and {V • [{jo5n + 2jono5j) < 
UtiUtj >]}\y = dxlijo^Ti ^ 2jorio5j) <^UtxUty >]+ < 
^ty > dy{j^5n^2jono5j)^{j^5n^2jono5j)dy{< u^y > 
)/2 =< u^y > [^ldy{5n)^2^QnQdy{5^)] since < Uty >= 
and dy{< u^y >) = 0, and {V • [(70 + 27ono(57) < 
UtiUtj >]}\z = since dz = and < Utz >= 0. The sec- 
ond term on the right side of Eq. [2] can also be omitted 
by taking the laser cycle average for po/jo. 

Then one can get: 
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And the last term F is not related to laser pressure and 
thermal pressure. 

Although the explicit relationship between a5a and 5n 
has not been found yet, from Eq. [2] one can see the 
thermal pressure can weaken the density modulation due 
to radiation pressure. In the region where the density 
is lower and the laser intensity is higher, the light pres- 
sure tends to expel the electrons in further, however, the 
thermal pressure tends to weaken such kind of plasma ex- 
pelling. With a proper transverse temperature, the bal- 
ance of thermal pressure and the perturbed light pressure 



can be built up just as shown in Figl3l Thus the excita- 
tion of the instability can be effectively suppressed. Since 
the thermal pressure is used only to balance the pertur- 
bation of the laser pressure which is small, a transverse 
temperature of MeV level is enough. We do not need the 
GeV level of temperature to compete with the total laser 
pressure. 

One may argue that the transverse temperature also 
causes transverse expansion of the target region. It seems 
the suppression of the instability would inevitably result 
into the target mass loss in the transverse direction. This 
is true. However, this kind of loss is far slower than 
the loss due to the instability. Moreover, recently Bu- 
lanov et al. have argued that this kind of transverse loss 
could actually benefit the remaining target acceleration, 
and results into the accelerated ions getting phase locked 
with the electromagnetic wave. This could, in principle, 
produce unlimited ion energy gain [12]. Certainly, the 
stable acceleration structure should be kept intact before 
the phase locking could actually occur. 

In summary, we have demonstrated that the proper 
transverse temperature can suppress the transverse insta- 
bility during the radiation pressure acceleration regime. 
We have shown that the surface erosion of the target can 
help in increasing the transverse electron temperature. 
This scheme keeps the target intact for a longer dura- 
tion compared with a normal flat target, resulting into 
the improvement of both final peak energy and spectrum 
quality of the accelerated ions. 

This work is supported by the DFG through TR-18 
project. MC acknowledges support by the Alexander von 
Humboldt Foundation and help from C.B. Schroeder on 
the relativistic warm plasma model. 



[1] T. Ditmire, et al, Nature 386, 54 (1997); J. Denavit, 
Phys. Rev. Lett. 69, 3052 (1992); A. Pukhov, ibid. 86, 
3562 (2001); B.M. Hegelich, et al, ibid. 89, 085002 
(2002); T. Esirkepov, et al, ibid. 89, 175003 (2002); T. 
Esirkepov, et al, ibid. 92, 175003 (2004); L. O. Silva, et 
al, ibid. 92, 015002 (2004); J. Fuchs, et al. Nature Phys. 
2, 48 (2006); Y. Lin, et al, Phys. Plasmas, 14, 056706 
(2007); M. Chen, et al, ibid. 14, 113106 (2007). 

[2] B.M. Hegelich, et al. Nature (London) 439, 441 
(2006) ;T. Toncian et al Science 312, 410 (2006). 

[3] M. Borghesi, et al, Phys. Plasmas 9, 2214 (2002). 

[4] E. Fourkal, et al, Med. Phys. 34, 577 (2007). 

[5] N. Naumova, et al, Phys. Rev. Lett. 102, 025002 (2009); 
M. Roth et al, Phys. Rev. Lett. 86, 436 (2001); 

[6] A. Macchi, et al, Phys. Rev. Lett. 94, 165003 (2005); X. 
Zhang, et al, Phys. Plasmas 14, 123108 (2007); A. P. 



[*] Present address: Lawrence Berkeley National Laboratory, 1 Cy- 
clotron Road, Berkeley CA, 94720. Email: MinChen@lbl.gov 



5 



L. Robinson, et a/., New J. Phys. 10, 013021 (2008); O. 

Klimo, et a/., Phys. Rev. Special Topics - Accelerators 

and Beams 11, 031301 (2008); X.Q. Yan, et a/., Phys. 

Rev. Lett. 100, 135003 (2008); A. Henig, et al, Phys. 

Rev. Lett. 103, 245003 (2009); A. Macchi, et a/., New J. 

Phys. 12 045103 (2010). 
[7] M. Chen, et a/., Phys. Rev. Lett. 103, 024801 (2009); M. 

Chen, et a/.. New J. Phys. 12, 045004 (2010). 
[8] F. Pegoraro, et a/., Phys. Rev. Lett. 99, 065002 (2007). 



[9] M. Chen, et a/., Phys. Plasmas 15, 113103 (2008). 
[10] A. Pukhov, J. Plasma Phys. 61, 425 (1999). 
[11] C. B. Schroeder, E. Esarey, Phys. Rev. E 81, 056403 
(2010). 

[12] S. V. Bulanov, E. Yu. Echkina, T. Zh. Esirkepov, L N. 
Inovenkov, M. Kando, F. Pegoraro, and G. Korn, Phys. 
Rev. Lett. 104, 135003 (2010). 



